Impurity Scattering Effect on Superconductivity and the Violation of Anderson 
Theorem in Ba(Fei ^.Co2)2As2 Single Crystals 
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Low-temperature specific heat (SH) and resistivity were measured on Ba(Fei_a;Coi:)2As2 single 
crystals in wide doping region. A sizeable residual specific heat coefficient 70 was observed in the 
low temperature limit of all samples. The specific heat jump near Tc, i.e. AC/T|tc and the upper 
critical field H^2 (H||c) were also determined. It is found that —70, AC/TjTc Hc2 ^^'^ 
shared a similar evolution with doping. All these can be well understood within the model of 
pairing symmetry when accounting the Co-dopants as unitary scattering centers in the FeAs 
planes. Our results give a direct evidence for the violation of the Anderson theorem in FeAs-based 
superconductors. 



Theoretically it was suggested that the superconducting 
pairing may be established via exchanging spin fluctua- 
tions between the electrons in the hole pockets (around F 
point) and the electron pockets (around M point), 3, i3| 
thus a model concerning s-wave symmetry with opposite 
sign between different bands (the so-called S^) was pro- 
posed. A direct evidence to prove this unique pairing 
manner is still lacking although some indirect evidence 
does indicate that the superconductivity vanishes grad- 
ually when the condition for this interpocket scattering 
deteriorates. 3 12 1 

In a conventional superconductor, the non-magnetic 
impurity will not lead to apparent pair-breaking effect, 
therefore no quasiparticle density of states (DOS) can 
be generated at the Fermi energy Yip. This was called 
as the Anderson theorem. [l^ In sharp contrast, in a d- 
wave superconductor, non-magnetic impurities can in- 
duce a high DOS due to the existence of nodes. Thus 
it is not strange when a large residual specific heat coef- 
ficient 7o (which measures actually the DOS at Ej?) was 
observed in the cuprate superconductors. [l9j| As for the 
case with the pairing symmetry of S^, it has been pointed 
out that non-magnetic irnpurity disorder could severely 
suppress Tc and the gap.[20| Recently, the effect of im- 
purity scattering in the case of S* pairing symmetry was 
considered and it was found that the DOS spectrum p{lo) 
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PACS numbers: 74.20.Rp, 74.70.Dd, 74.62.Dh, 65.40.Ba 

The discovery of superconductivity in iron pnictides 
have generated enormous interests in the community of 
condensed matter physics. fH One of the key issues here 
is about the superconductivity mechanism. Phcnomcno- 
logically it has been found that in many, if not all, struc- 
tures of the iron-pnictidc materials, the parent phase has 
a long range antiferromagnetic (AF) ordered state, and 
the superconductivity is induced by suppressing this AF 
order.[^|3,[^, Q Another important issue concerning the 
superconductivity mechanism is about the pairing sym- 
metry. Up to date, experimental results gave rather con- 
tradicting conclusions about the pairing symmetry in the 



can be significantly modified leading to a finite DOS at 
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22, |23|, |2J] Conclusions about the residual DOS 
from thermal transport measurement are controversial 
with each other. [H, [H, [l^, ,28] To verify this theoreti- 



iron pnictide superconductors. 0, S i 0, [nl. 111, Q AC/T|t, and 



cal hypothesis, specific heat (SH) measurement may be 
a good choice because it is straightforward to get the in- 
formation of DOS at the Fermi level. In this Letter, we 
report the low-temperature SH under different magnetic 
fields on the Ba(Fei_2;Coa;)2As2 single crystals from un- 
derdoped to overdoped region. We found a sizeable value 
of 7o for all samples in the low temperature limit. The 
doping dependence of 70 anti-correlates with that of T^, 
1%. These behaviors were well inter- 



preted within the model of S pairing symmetry. 

The Co-doped Ba(Fei_2,Coa;WAs2 single crystals were 
grown by the self-fiux method[17|. The samples for the 
SH measurement have typical dimensions of 2.5 x 1.5 
X 0.2 mm^. The dc magnetization measurements were 
done with a superconducting quantum interference de- 
vice (Quantum Design, SQUID, MPMS7). The resistiv- 
ity and the specific heat were measured with a Quantum 
Design instrument physical property measurement sys- 
tem (PPMS) with the temperature down to 1.8 K and 
the magnetic field up to 9 T. We employed the thermal 
relaxation technique to perform the specific heat mea- 
surements. To improve the resolution, we used a latest 
developed SH measuring frame from Quantum Design, 
which has negligible field dependence of the sensor of the 
thermometer on the chip as well as the thermal conduc- 
tance of the thermal linking wires. 

In Fig. 1(a), we show the temperature dependence of 
resistivity under zero field for five samples with nominal 
doping compositions of a; = 0.06, 0.07, 0.08, 0.12, and 
0.15, respectively. The sample with x = 0.08 was found 
to be optimally doped with the highest onset transition 
temperature T°"'*'^* « 24.5 K. In the underdoped region 
{x < 0.08), an upturn in the resistivity curve above Tc 
can be easily seen, which was supposed to be related with 
the structural and antiferromagnetic ( AF) transition. 
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FIG. 1: (color online) (a) Temperature dependence of resistiv- 
ity for the Ba(Fei-xCoa;)2As2 single crystals in wide doping 
range under zero field. (b)The enlarged view of the resistivity 
data near the superconducting transition. (c)The dc mag- 
netization data measured with H = 20 Oe for the zero field 
cooling (ZFC) process. The curves were normalized by the 
magnetization data obtained at 2 K. 



An enlarged view of the p{T) curves near Tc was shown 
in Fig. 1(b). We also measured the dc magnetization 
of the samples, which was displayed in Fig. 1(c). The 
rather sharp transitions suggest the high quality of our 
samples. 

We show the raw data of SH for the sample x = 
0.08 in the main frame of Fig. 2. The red solid line 
displays roughly the tendency of the normal state SH, 
Cnorm/T = 7„ -l- Cph/T, where 7„ is the electronic con- 
tribution and Cph/T is the phonon contribution obtained 
based on a simple polynomial fit in the normal state. 
This will not be relied on to analyze our data. A clear 
anomaly due to the superconducting transition can be 
observed at about 23 K in the zero-field data. A mag- 
netic field of 9 T suppresses the anomaly remarkably and 
also moves the transition to lower temperatures. We also 
show the enlarged view of the data in the low-T region 
in Fig. 2(a). One can see the roughly linear behavior 
in the C/T vs plot in low-T region. Surprisingly, no 
clear Schottky anomaly was detected in all the samples, 
which may suggest that the Co-doping here induces no 
local magnetic moment which would on the other hand 
give a large contribution to SH as the Schottky anomaly. 
It is clear that the magnetic field enhances the low-T 
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FIG. 2: (color online) Main frame: The raw data of SH for 
the optimally doped sample {x = 0.08) under different fields. 
The red solid line is guide to eyes which displays the rough 
tendency of the normal state SH (see text). The inset (a) 
shows an enlarged view of the same data in the low-T region. 
Inset (b) shows the field dependence of the field-induced term 
A7(if) = [C{H) - C(0)]/T at K obtained by linearly ex- 
trapolating the low-T data to zero temperature. The green 
solid, and red dash-dotted lines are the fit to the clean d-wave 
prediction A'y{H) = A^/H, and the dirty d-wave prediction 
A-f{H) = A{H/Hc2)los[B{H,2/H)], respectively. 



SH continuously, indicating the increase of quasiparticle 
DOS at Fip induced by magnetic field. We will discuss 
this issue later. 

In order to have a comprehensive understanding, we 
also measured the temperature and field dependence of 
SH on samples in wide doping region. We extracted the 
SH difference between T and 9 T and showed the re- 
sults in Fig. 3(a). From the main frame of Fig. 2, we 
can see that a magnetic field of 9 T can not suppress 
the superconductivity completely, but it shifts the super- 
conducting transition to a distinguishable lower temper- 
ature. As a result, we can evaluate the height of the SH 
anomaly AC/T\t^ near from the difference of C/T at 
T and 9 T. It is clear that the optimal doped sample 
with highest has the largest anomaly AC/Tlr^ « 28.6 
mJ/mol K^. This value is quit e comparable with that 
reported by other groups. (29ll30l [sif In each doping side 
(underdoping or overdoping), AC/T\t^ seems to display 
a monotonic increase with Tc- This behavior is quali- 
tatively consistent with that reported in Ref. 3]| where 
a scaling behavior of AC/T\t^ oc (Tc)^ was reported. 
However, we note that there is a clear difference between 
the underdoped and the overdoped regions. For example, 
the sample with x = 0.07 has a higher Tc while showing a 
smaller AC /T\t^ compared with that of the sample with 
X = 0.12. As will be discussed below, we attribute this 
difference in the underdoping and overdoping regions to 
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FIG. 3: (color online) (a) The difference of SH between 
T and 9 T for samples in wide doping region, (b) Low- 
temperature SH data for samples with different doping levels 
under zero temperature. One can see a nonmonotonic evolve- 
ment of the residual value 70 with doping. The departure 
from linear behavior at about 7 K for the sample with x — 
0.15 was caused by the superconducting transition. 



the different mechanism that governs the evolvement of 
AC/T\t, with Tc. 

In Fig. 3(b), we present the low-T SH data at zero 
field. A linear extrapolation of the low-T data finds im- 
mediately that there is a sizeable value of the residual SH 
coefficient 70 for all samples. A closer scrutiny realizes 
that 7o has a nonmonotonic doping dependence on the 
Co-doping concentration x. We must stress that the size- 
able value of 7o found in present samples should not be 
simply attributed to the non-superconducting fraction. 
The reasons are as following: (1) A minimum of 70 was 
observed just at the optimal doping point (see below). 
From the chemistry point of view, however, there is no 
reason to believe that the non-superconducting fraction 
should be the lowest in the optimally doped sample; (2) 
Counting the magnetization signal in the low tempera- 
ture region finds that the magnetic shielding is beyond 
95 % for all the samples. 

In order to clarify the origin of 70 in present samples, 
we have extrapolated the low-T SH data shown in inset 
(a) of Fig. 2 to zero temperature linearly and obtained 



the field-induced term A7(ff) = [C{H) - C(0)]/T at 
K, which was shown in inset (b) of Fig. 2 for the opti- 
mally doped sample. One can see that Aj{H) rises up 
quickly and shows a nonlinear tendency below about 1 
T. Whereas it displays a clear linear behavior above 1 
T. Similar feature was observed in samples with other 
doping levels. We at first attempted to fit the data with 
the relation Aj{H) — A^/H predicted for d-wave sym- 
metry in the clean limit. [l^] The result was displayed 
by the green solid line in the inset (b) of Fig. 2. It is 
clear that this fitting curve cannot describe the exper- 
imental data at all. Secondly, we fitted our data us- 
ing the relation for d-wave superconductors in the dirty 
limit, Aj{H) = A{H/H^^)log[B{H^^/H)].^] Here B 
is a constant which approximates 7.26 for a triangular 
vortex lattice. We left A and ^-s the free fitting pa- 
rameters. The best fitting result was shown by the red 
dash-dotted line. Again this curve departs from the ex- 
perimental data, especially it cannot reflect the kink fea- 
ture around 1 T and the linear feature above 1 T. So 
we can exclude the presence of the superconducting gap 
with d-wave symmetry, either in the clean or dirty limit. 
Consequently, the finite DOS found in the present system 
cannot be attributed to the impurity scattering effect for 
a d-wave superconductor. 

The doping dependence of the extracted 70, along with 
Tc, was shown in Fig. 4(a). The curve of Tc vs x 
formed an asymmetric dome, while the 70 vs x curve 
showed an anti-correlated behavior. We argue that this 
behavior can be explained by the Co-induced impurity- 
scattering effect. Numerical calculations using the T- 
matrix method have shown that, in a superconductor 
with S* pairing symmetry, the fully open ed g ap of a clean 
state will be filled up by impurity states. 2lJ Therefore a 
finite DOS at Ef may rise up (forming the so-called gap- 
less state) if the scattering strength becomes stronger. In 
the unitary limit, the residual SH coefficient 70 may be 
expressed by the impurity concentration riimp and the su- 
perconducti ng g ap in a simple form 70 cx [riimp/ A^)", 
with Of > 0. 3J] The index a approximates 0.5 for a d- 
wave superconductor, [ssj 

In the underdoped region, assuming a proportionality 
between riimp and the number of Co-dopant, since the 
magnitude of A^ increases with x more rapidly (see dis- 
cussion later on Fig. 4(d)) than a linear increase of Uimp, 
therefore 70 was reduced with the increase of doping, 
one thus qualitatively understands that 70 anti-correlates 
with As c>c Tc- Actually the realistic case is more com- 
plicated: the superconductivity and AF states compete 
with each other in the underdoped region, and there will 
be less and less contributions of DOS given by the AF 
state in the T = K approach with adding more Co- 
dopants into the system. These two factors lead to the 
dropping-down behavior of 70 versus x. 

In the overdoped region, the AF order was suppressed 
completely. However, as we have addressed, the Co- 
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FIG. 4: (color online) Doping dependence of (a) the resid- 
ual SH term 70, (b) minus of the residual term —70, (c) the 
superconducting SH anomaly AC/TItc and (d) square root 



of the upper critical field 
comparison. 



plotted along with Tc for 



doping will deteriorate the spin-fluctuations and weaken 
the pairing strength, resulting in the decrease of A 
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Meanwhile 



''imp 



(oc x) keeps rising. These two factors 



lead to the quick increase of 70 with doping in the 
overdoped region. In Fig. 4(b) we also showed —70 to- 
gether with Tc. Surprisingly one can see a quite good 
consistency between the doping dependence of —70 and 
Tc- This good consistency is understandable because 70 
reflects how many spin-fluctuation-mediated scattering 
channels, which are responsible for the Cooper pairing, 
are blocked away by the impurities. 

As for the doping dependence of the SH jump 
AC/TIt^, which were shown in Fig. 4(c), we can ex- 
plain it based on the variation of 70 with x. The BCS 
theory tells that the height of SH jump is proportional 
to the effective normal state SH coefficient 7e/ / : 



AC 



It, oc 7e// oc (7; 



bare 



(1) 



Here 7^°'"'^ is the bare value of SH coefflcient in the normal 
state (weakly dependent on doping [35|]) if the impurity 
scattering and the competition of the AF state would not 
exist, and /3 > 1. The impurity-scattering may modify 
the DOS spectrum p{lo) in the w A approach and 
suppress the height of the SH jump AC/T\t^, giving a 
value of P larger than 1. From equation (2), we can 
expect a roughly consistent tendency between AC/T\t^ 
and —70. Recalling the fact that Tc correlates linearly 
with —70, one can easily see the reason for the similar 
evolvement tendency between Tc and AC/T\t^ with x, 
as shown in Fig. 4(c). 

In above discussion, we have shown that the pa- 
rameters —70, AC/T\t^ and Tc share a similar dop- 



ing dependence. The key player here is actually As 
which is estimated in the following. In Fig. 4(d), 
we show the doping dependence of the square root of 



the upper critical fleld (H||c-axis) 



along with 



Tc, where Hc2 was determined using the Werthamer- 
Helfand-Hohenberg relation [3^ from the field dependent 
resistivity data (not shown here) . One can see that both 
set of data overlap quite well. This is understandable be- 
cause the Ginzburg-Landau theory has given the relation 



$0 



^TT^c^ab 



OC 



(2) 



where $0 is the flux quantum and ^c (Cab) is the coherence 
length in the direction of c axis (ab plane). As a result, 
y^H^2 is proportional to Ag oc Tg. 

In summary, we studied the low-temperature SH and 
resistivity on Ba(Fei_a;Coa;)2As2 single crystals in wide 
doping region. A sizeable residual SH coefficient 70 in 
the low-T limit and clear SH jumps were detected in all 
samples. It is found that —70, AC/T\t^, y^S^, and 
Tc all share a similar evolution with doping amount x. 
All these behaviors were interpreted within the model of 
S^ pairing symmetry considering the Co-doping induced 
scattering effect in this system. 
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